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Abstract 

Dynamics of attractive ultra-cold bosonic clouds in one dimension are studied by solving the 
many-particle time-dependent Schrodinger equation. The initially coherent wave-packet can dy- 
namically dissociate into two parts when its energy exceeds a threshold value. Noticeably, the 
time-dependent Gross-Pitaevskii theory applied to the same initial state does not show up the 
splitting. We call the split object fragmenton. It possesses remarkable properties: (1) it is two-fold 
fragmented, i.e., not coherent; (2) it is dynamically stable, i.e., it propagates almost without dis- 
persion; (3) it is delocalized, i.e., the two dissociated parts still communicate with one another. A 
simple static model predicts the existence of fragmented states which are responsible for formation 
and dynamics of fragmentons. 

PACS numbers: 03.75.Kk, 05.30.Jp, 03.75.Nt, 03.65.-w 
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Dynamics of attractive low- dimensional dilute Bose gases have attracted much attention 
. The formal equivalence of the Gross-Pitaevskii (GP) equation used to 
describe these quantum gases and the non-linear Schrodinger equation used in non-linear 
optics has led to transfer of knowhow and, most importantly, to the prediction of bright 
solitons in attractive Bose gases 

US. 

More recently, fascinating and challenging exper- 
iments on bright matter-wave solitons have been reported [3, 4| and explained using GP 
theory, making attractive Bose gases a cornerstone of GP theory. Nowadays, the GP theory 
is widely utilized to study attractive dilute Bose gases. 

GP theory is a mean- field approximation to the quantum many-boson problem, in which 
the Bose system is assumed to remain condensed during its evolution in time. Recently, the 
general validity of this assumption has been questioned in the theoretical study on dynamical 
depletion of attractive one-dimensional (ID) Bose systems at non-zero temperatures 
In the present Letter we consider the dynamics of attractive ID Bose systems from the 
full many-body point of view and relate the many-body and GP results. We shall see 
that the dynamics of attractive ID Bose gases can be much richer than anticipated before 
because along with condensed states the system supports a new kind of low-lying excited 
states which are fragmented, and thereby not available within the framework of GP theory. 
These low-lying fragmented excited states are responsible for the formation and evolution 
of dynamically-stable fragmented objects in attractive ID Bose systems. These objects are 
distinct from the solitons discussed above. 

We solve the time-dependent many-body Schrodinger equation for N attractive bosons 
in ID interacting via contact potential U(x — x') = \o5(x — x') in free space within 
the framework of the recently developed multi-configurational time-dependent Hartree 
for bosons (MCTDHB) jjj. For comparison we also integrate the respective time- 
dependent GP equation. In GP theory the wavefunction has the form *$?(xi, . . . ,x^,t) = 
4>i(xi, t)<pi(x2, t) ■ ■ ■ 4>i(xn, t), i.e., all bosons reside in a single orbital 4>i(x, t). MCTDHB(M) 
goes much beyond GP in that the bosons are distributed over M orbitals {4>i(x, t)} and all 
possible distributions are considered. The MCTDHB(M) ansatz reads ^(xx, X2, ■ ■ ■ , xn, t) = 

EmCmW^l^l-i) •••<h( 

t) ■ ■ ■ 4> m {xn, t) where S is symmetriza- 
tion operator and m = (m^mj, ■ ■ - tum) occupation numbers. Within this theory a quanti- 
tative description of the time evolution of bosonic systems is achieved by optimizing all M 
orbitals used to construct the many-body expansion and the expansion coefficients them- 
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selves at each point in time utilizing time-dependent variational principle j9(. 

At t=0 we take as initial condition a totally condensed atomic cloud of a very-well defined 
shape and propagate the many-body MCTDHB(2) and GP equations in time and compare 
the obtained densities. In this work we take initial coherent wave-packets where all bosons 
are in either Sech [jx]- or Gaussian-shaped normalized functions. These choices represent 
different ways to generate the initial conditions, leading, however, to similar physical results. 
Varying the widths of the initial wave packet and keeping the number of bosons and their 
interparticle interaction strength fixed allows us to change the energy of the initial wave 
packet and to activate thereby different excited states of the many-boson attractive system. 
In this study we allow only for initial wave-packets which energetically do not exceed the 
formation energy of two-hump soliton trains. Thus, by construction all the results observed 
in the present study cannot be attributed to soliton trains. 

In Fig. [1] we present the results of the many-body Schrodinger equation for N=1000 
attractive bosons with Ao = —0.008 for several 5"ec/i-shaped initial wave-packets of different 
widths. Here and in the following we work in dimensionless units which are readily arrived 
at by introducing a convenient length scale L (say, the scattering length) and dividing the 
Hamiltonian by -^rs, where m is the mass of a boson. The densities are plotted as a function 
of time. The many-body MCTDHB densities depicted in the left column can be compared 
with the respective time-dependent GP densities plotted in the right column obtained for 
the same system with identical initial conditions. The ground eigenstate of the many-boson 
attractive system in free space at the GP level of description is known analytically - it is the 
famous Sech [jx] with optimal exponent 7 = l*°l^ r ~ 1 ) . For the system under investigation 
the optimal exponent corresponding to the ground eigenstate is 7 ~ 4.0. In the first study we 
take a slightly broader initial 5ec/i-function with 7/ = 3.0. The corresponding many-body 
and GP dynamics are depicted in the upper left and right panels of Fig. [TJ respectively. 
We see that the width of the wave-packet changes in time rather periodically describing 
thereby "breathing" of the attractive atomic cloud. The many-body and GP dynamics are 
quite similar (but not identical) confirming the solitonic character of the ground state also 
at the many-body level of description. In the second study presented in the middle panels 
of Fig. [T]we compare many-body and GP dynamics for a broader initial <Sec/i-function with 
77/ = 1.2. We see that dynamics of the first five "breathing" oscillations are quite similar 
within both theories, while at longer times the many-body theory shows attempts of the 
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density to split into two parts. The bottom panels of Fig. [T] show the results of the third 
study, where we compare many-body and GP propagations of an even broader initial Sech- 
function with 7/// = 1.0. Again, the first three "breathing" oscillations of the GP dynamics 
resemble the respective many-body ones, then, however, the many-body dynamics reveal 
a fascinating feature of the evolution - the initially condensed single Sech- shaped cloud 
dissociates into two smaller equal clouds which move symmetrically apart from each other 
with approximately constant velocity. 

We stress several points. First, the dynamical many-body splitting is a general phe- 
nomenon found for many initial conditions and for wave-packets of other shapes as well, 
e.g., for symmetric and asymmetric Gaussian-shaped profiles. Asymmetric initial packets 
dissociate into unequal parts. Second, the phenomenon takes place only if the energy of the 
initial cloud exceeds some threshold value (see below). Third, by increasing the widths of the 
initial wave packets we increase their total energies making thereby different excited states 
of the quantum many-boson system energetically accessible. The existence of the splitting 
phenomenon found only at the many-body level of description indicates that in attractive 
ID Bose systems there is a new kind of low-lying excited states, which is not available within 
the framework of GP theory. 

To get a deeper insight into the physics, we investigate next the many-body structures of 
the evolving wave-packets in some details. By diagonalizing the reduced one-body density 
matrix at each time step we observe how its eigenvalues, i.e., natural occupation numbers 
(rii) and its eigenvectors, i.e., natural orbitals evolve during the dynamics. The evolutions 
of the natural occupation numbers for all the three initial conditions reported in Fig. [1] are 
plotted in Fig. [2j We choose log scale to plot the time evolution of the Ui/N in %. At t=0 all 
the bosons of the initial wave-packets reside in one and the same natural orbital indicating 
that these states, according to the usual definition [lOj are condensed. In time we see that 
another orbital acquires some population revealing the depletion of the initially-condensed 
state. Interestingly, during the times where the system shows "breathing" dynamics, the 
changes of the respective natural occupations are moderate. Indeed, for 7/ = 3.0 where 
GP theory predicts very similar breathing dynamics, the occupation of the second natural 
orbital does not exceed 2%. In the other studies (77/ = 1.2,7/// — 1-0) the second nat- 
ural occupation number oscillatory increases until 10-15% as long as the density exhibits 
"breathing" oscillations. Then, the system fragments (see Fig. [2]). 
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Let us concentrate on the third study. As time proceeds, the occupation of the second 
natural orbital increases and around the point where the cloud starts to dissociate into two 
parts, the two natural occupation numbers become equal and macroscopic. According to the 
usual definition such a system is two-fold fragmented. Clearly, GP theory is inapplicable 
here. Further evolution reveals oscillations of the occupation numbers around this ideally 
two-fold fragmented state. This results from the fact that the respective natural orbitals 
themselves are delocalized objects, so that even at very large separations the dissociated 
parts remain connected and can "talk" to each other (see below). The two separated packets 
propagate almost without dispersion like ordinary solitons (see also Fig. [1]). We call this 
new physical object fragmenton because it combines macroscopic fragmentation of the wave 
function and dynamical properties of a soliton. We recall that in contrast to the ordinary 
soliton or soliton train solutions which are coherent objects, fragmenton is fragmented, i.e., 
characterized by macroscopic occupations of several (here two) natural orbitals. 

For the sake of interpretation, we construct a simple static model which is able to describe 
two-fold fragmented states \ni,n 2 ), where n\ bosons reside in 0i and n 2 = N — n\ in <p 2 
orbitals. As the orbitals constituting this state we choose delocalized gerade and ungerade 
superpositions of two identical Sec/i-functions with equal exponent 7 m placed at ±A with 
respect to the origin: 

0i,2 = *i,2 [Sech(j m (x - X )) ± Sech{-f m (x + X ))) , (1) 

where Yi j2 are the respective normalization factors. Our goal is to find within this ansatz the 
lowest energy eigenstates of the N identical interacting bosons. We consider the exponent 
7 m and position Xq of the S'ec/i-functions as variational parameters. Additionally we can 
change the occupation number n\ (n 2 = N — ni). The expectation energy of a two-fold 



fragmented state is known 121 ]: 



E{rii) = riihu + A — — l - / \4>i\ 4 dx + n 2 h 22 



+x Mn 2 1) ^ dx + 2X onin2 I Ifaflfafdx (2) 



where ha = {<pi\ — \-£s\<j>i) are the expectation values of the kinetic energy operator. The 
minima of this energy functional are plotted in Fig. [3]for our system of N=1000 bosons with 
Aq = —0.008 as a function of n\. 
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The ground state corresponds to the situation where all the bosons are condensed in one 
localized orbital \N, 0). For an open system this is exactly the bright soliton GP solution 
with X = and 7 m = l^°K^zM ) which is reproduced by our model. In the right lower 
corner of Fig. [3] we schematically plot this function. If we increase the population of the 
second orbital, the total energy of the system increases and the density starts to split, i.e., 
X has non-zero values. For example for rii = 3N/4, the minimum energy E/N = —1.7135 
is obtained for individual 5ec/i-functions slightly (X = 0.2907) separated from each other 
and have exponent 7 m = 2.9839. In Fig. [3] we depict this energy and respective orbitals 
\Jriil N(f)i . However, our simple model accounts for two branches of physically different 
solutions. For any given ri\ - there is another branch of solutions which is much higher 
in energy and characterized by large values of optimal X , i.e., by two very-well separated 
humps. Indeed, for n\ = 3N/4 such a solution has optimal energy E/N = —1.2586 with 
exponent 7 m = 2.7480 and separation X = 4.0. We depict the respective energy point and 
pair of orbitals in Fig. [3j Reducing rii further, we increase the energy of the lower branch and 
decrease the energy of the upper branch which approaches its minimum at rt\ = n 2 = N/2. 
At this point, as we can see in Fig. [3], one encounters a bifurcation of the branches. Increasing 
the occupation of the second orbital further, the energy increases until at n\ = 0, n 2 = N 
we arrive at the two-soliton anti-phase (ungerade) solution, plotted in the left upper corner 
of Fig. [31 

There are three points of relevance indicated in Fig. [3j Ground state (GS), excited 
state (ES) and the bifurcation point (BF) which corresponds to the ideally fragmented 
(rii = ri2 = N/2) state of minimal energy. In our model the optimal energies of all these 
states have a unique form E state = — Nj state 2 /6, with different optimal exponents Estate- 
Indeed, the ground state energy of the bright soliton is known: Eq P = —N'-fcs 2 /6 with 
7G5 = | Ao| (iV — l)/2. The energies of the coherent two-soliton gerade and ungerade states 
can be obtained by minimizing the GP energy functional assuming zero overlap of the 
constituting solitons. Then these degenerate energies are Eqp = —N^es 2 /& with exponent 
Ies — W\{N — l)/4. Finally, assuming that the Sech-i unctions forming the orbitals of 
ideally fragmented state do not overlap, we get by minimizing Eq. ([2]) the optimal exponent 
Ibf = |Ao|(3iV/2 — l)/4 and bifurcation energy E BF = —N~f BF 2 /6. From these energies 
we see that the ideally two-fold fragmented state is separated from the ground state by 
the energy gap AEfr w 7iV 3 Ao 2 /384. This energy can be viewed as the threshold for the 
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activation of the fragmenton. The energy difference between the ground and excited GP 
states is AEes ~ 12iV 3 Ao 2 /384. Therefore, in attractive one-dimensional Bose gases two- 
fold fragmented delocalized states are always energetically more favorable than the respective 
two-soliton coherent gerade and ungerade states. 

Now we can interpret the wave-packet evolutions presented in Fig. [T] in terms of the states 
of the system (see Fig. [3]). If the energy of the initial wave-packet relative to the ground 
state does not exceed the fragmenton activation energy AEfr, the model permits dynamics 
only within the states of the lowest branch. Indeed, in the first study the energy Ej of 
the initial wave packet (see Fig. [3]) is insufficient for activating the fragmenton. The many- 
body scenario plotted in the left upper panel of Fig. [T] shows small-amplitude "breathing" 
oscillations around the ground state. In contrast, if the activation threshold is overcome, as 
it happens in the second (Ejj) and in particular in the third (Em) study, the fragmenton 
channel becomes accessible energetically and strongly affects the dynamics. This is reflected 
in the left middle and lower panels of Fig. [T] where the initial wave packet attempts to split 
and successfully dissociates into two parts, respectively. The respective Ej, En and Em 
energies are shown in Fig. [3j 

Let us summarize. We explore the challenging many-body dynamics of attractive ultra- 
cold bosonic clouds in ID by varying the width and thereby the total energy of the initial 
coherent S'ec/i-shaped wave-packet. We solve the many-boson time- dependent Schrodinger 
equation and show that the many-body dynamics of attractive condensates can differ con- 
siderably from the respective GP dynamics. We find that when the initial energy of the 
condensate exceeds some threshold, the condensed cloud at the many-body level of descrip- 
tion tends to split and can dynamically dissociate into two separate parts, while at the 
GP level it reveals only "breathing" dynamics and remains unsplit. To understand this 
phenomenon we first analyze the many-body wave-function and attribute the split clouds 
to a two-fold fragmented delocalized state. We call this new class of dynamical solutions 
fragmentons because on the one hand they propagate almost without dispersion, and on the 
other hand, in contrast to solitons, they are not coherent objects. A simple static model is 
introduced to show that the two-fold fragmented state constructed from in- and anti- phase 
two-soliton-like functions is energetically much more favorable than coherent two-soliton so- 
lutions. To estimate the window for fragmenton formation we provide simple estimations of 
the energies of the involved states in terms of the number of particles iV and their interaction 
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strength X . 

Our findings imply that the dynamics of attractive Bose gases in low dimensions is much 
richer than anticipated before and we hope that our work will stimulate experiments on 
fragmentons. 
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FIG. 1: (Color online) MCTDHB many-body (left) vs. GP (right) dynamics for N=1000 attractive 
(Ao = —0.008) bosons in one-dimension. The densities are plotted as a function of time for initially- 
coherent wave-packets Sech[jx] of different widths. While the GP dynamics shows "breathing" 
oscillations for all initial widths, the many-body dynamics changes dramatically from "breath- 
ing" oscillations (7/ = 3.0), to attempts for dynamical splitting (777 = 1.2), and to dynamical 
dissociation into two parts (7/// = 1.0). All quantities are dimensionless. 
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FIG. 2: (Color online) Evolutions of natural occupation numbers n\/N and U2/N of the reduced 
one-body densities plotted in the left panels of Fig. [Jon a log scale in %. At t=0 all initial states 
are condensed: ri\/N = 100%. In the first study (7/ = 3.0) the dynamics describe "breathing" 
oscillations only. In the second (7/7 = 1.2) and third (77// = 1.0) studies the systems evolve to be 
two-fold fragmented after few "breathing" oscillations. Notice that in GP n\/N is always 100 %. 
All quantities are dimensionless. 
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FIG. 3: (Color online) Energy diagram E(n\) of the model system introduced in Eqs. (jl|2p for 
two-fold fragmented \n\ : ri2) states n\ + n<i = iV. There are two branches of solutions bifurcating 
at m = iV/2 (bifurcation energy Ebf)- Shown are also the orbitals w rn /N(f>i corresponding to 
n\ = 3N/A at both branches and to the GP solutions of n\ = N and ri2 = N at energies Eqp and 
Eq P - Ej, Ejj, Ejjj are the energies of the three cases studied in Fig. [TJ To activate a fragmenton 
the energy of the initial cloud must exceed Ebf- All quantities are dimensionless. 
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